Radiating black holes in Einstein- Yang-Mills theory and cosmic censorship 
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Exact nonstatic spherically symmetric black-hole solution of the higher dimensional Einstein- 
Yang-Mills equations for a null dust with Yang-Mills gauge charge are obtained by employing 
Wu-Yang ansatz, namely, HD-EYM Vaidya solution. It is interesting to note that gravitational 
contribution of YM gauge charge for this ansatz is indeed opposite (attractive rather than repul- 
sive) that of Maxwell charge. It turns out that the gravitational collapse of null dust with YM 
gauge charge admit strong curvature shell focusing naked singularities violating cosmic censorship. 
However, there is significant shrinkage of the initial data space for a naked singularity of the HD- 
Vaidya collapse due to presence of YM gauge charge. The effect of YM gauge charge on structure 
and location of the apparent and event horizons is also discussed. 

PACS numbers: 04.20.Jb, 04.70.Bw, 04.40.Nr 



I. INTRODUCTION 



It is rather well established that higher dimensions pro- 
vide a natural playground for the string theory and they 
are also required for its consistency. In fact, the first suc- 
cessful statistical counting of black-hole entropy in string 
theory was performed for a five dimensional black-hole 
[l|. This example provides the best laboratory for the 
microscopic string theory of black holes. The classical 
general relativity in more than four space-time dimen- 
sions has also been the subject of increasing attention 
0|. Even from the classical standpoint, it is interesting 
to study the higher dimensional (HD) extension of Ein- 
stein's theory, and in particular its black-hole solutions 
[l| . The physical properties of these solutions have been 
widely studied. Interest in HD black holes has been fur- 
ther intensified in recent years, thanks to the well known 
AdS/CFT correspondence which envisions a correspon- 
dence between string theory (or supergravity) on asymp- 
totically locally anti-de Sitter backgrounds and the large 
TV limit of certain conformal field theories defined on the 
boundary-at-infinity of these backgrounds . There 
seems to be a general belief that endowing general rel- 
ativity with a tunable parameter namely the space-time 
dimension should also lead to valuable insights into the 
nature of the theory, in particular into its most basic ob- 
jects: black holes. For instance, four dimensional (4D) 
black holes are known to have a number of remarkable 
features, such as uniqueness, spherical topology, dynam- 
ical stability, and the laws of black-hole mechanics. One 
would like to know which of these are peculiar to four 
dimensions, and which are true more generally? At the 
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very least, such probings into HD will lead to a deeper 
understanding of classical black holes and of what space- 
time can do at its most extreme. There is a growing 
realization that the physics of HD black holes can be 
markedly different, and much richer than in four dimen- 
sions 0- 

However, there are very few nonstatic black- hole solu- 
tions known; one of them is the Vaidya black-hole solu- 
tion. The Vaidya solution Q is a solution of Einstein's 
equations with spherical symmetry for a null dust (radi- 
ally propagating radiation) source (a Type II fluid) de- 
scribed by energy-momentum tensor, T f, = ipn a nb, n a 
being a null-vector field. Vaidya's radiating star met- 
ric is today commonly used for two purposes: (i) as a 
testing ground for various formulations of the Cosmic 
Censorship Conjecture (CCC) and (ii) as an exterior so- 
lution for models of objects consisting of heat-conducting 
matter. Recently, it has also been employed with good 
effect in the study of Hawking radiation, the process of 
black-hole evaporation and in the stochastic gravity 
program pjj. It has also advantage of allowing a study 
of the dynamical evolution of horizon associated with a 
radiating black- hole. Further, various situations involv- 
ing spherically symmetric source as a mixture of a perfect 
fluid and null dust have been studied both in 4D [H Q3] 
as well as in HD [13h19| . 

It is of interest to consider models based on different in- 
teracting fields including the Yang-Mills. In general, it is 
difficult to tackle Einstein- Yang-Mills (EYM) equations 
because of the nonlinearity both in the gauge fields as well 
as in the gravitational field. The solutions of the classi- 
cal Yang-Mills fields depend upon the particular ansatz 
one chooses. Wu and Yang [2Cj] found static spherically 
symmetric solutions of the Yang-Mills equations in flat 
space for the gauge group SO (3). A curved space-time 
generalization of these models has been investigated by 
several authors (see, e.g., [ll[). Indeed Yasskin [2l| has 
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presented an explicit procedure based on the Wu-Yang 
ansatz [2(| which gives the solution of EYM rather triv- 
ially. 

Using this procedure, Mazharimousavi and Halilsoy 
[22Tj2~3 | have found a sequence of static spherically sym- 
metric HD-EYM black-hole solutions. The remarkable 
feature of this ansatz is that the field has no contribu- 
tion from gradient; instead, it has pure YM non-Abelian 
component. It, therefore, has only the magnetic part. It 
would be interesting to study the effect of pure YM field 
on gravitational collapse of null dust. That is the main 
motivation of this paper. 

By employing the Wu-Yang ansatz, we shall present a 
class of HD nonstatic solutions describing the exterior of 
radiating black holes with null dust endowed with gauge 
charge. That is, we find analogue of the HD-Vaidya solu- 
tion in EYM theory. We shall also consider the effect of 
YM gauge charge on the collapse of null dust described by 
the Vaidya solution onto a flat Minkowski cavity in HD. 
In the next section, we write the effective EYM equation 
in HD and find the generalized Vaidya solution namely 
HD-EYM Vaidya solution which would be followed in 
Sec. Ill by discussion of energy conditions and horizons. 
In Sec. IV, we shall study the gravitational collapse of 
null dust with YM gauge charge and then we conclude in 
Sec V. 



II. HD-VAIDYA LIKE METRIC IN 
EINSTEIN- YANG-MILLS THEORY 

We consider (N — 1)(N — 2)/2 parameter Lie group 



with structure constant C, 



(«) 

why 



The gauge potentials 



A { a a) and the Yang-Mills fields are related through 
the equation 



= d a A 



(a) 
b 



Then one can choose the gravity and gauge field action 
(Einstein- Yang-Mills), which in TV-dimensions reads: 



dx 



N 



M 



R 



(JV-l)(JV-2)/2 

\ " p{a)p(a)ab 



a = l 



ab 



(2) 

Here, g = det(g a b) is the determinant of the metric ten- 
sor, R is the Ricci Scalar [22LI231]. and A„ are the gauge 
potentials, ijv is the action of null dust. We note that 
the internal indices {a, /3, 7, ...} do not differ whether in 
covariant or contravariant form. We introduce the Wu- 
Yang ansatz in iV-dimension [221-12"? 



as 



Q 



(yX^diXj x j dx 



(3) 



2 < i<N-l, 
1 < j < i ~ 1, 

!<(«)< (iV-l)(iV-2)/2, 



where the super indices a is chosen according to the val- 
ues of i and j in order and we choose a = Q [22T4241 ] . 
The Wu-Yang solution appears highly nonlinear because 
of mixing between space-time indices and gauge group 
indices. However, it is linear as expressed in the nonlin- 
ear gauge fields because purely magnetic gauge charge is 
chosen along with position dependent gauge field trans- 
formation [2l|. The YM field 2 form is defined by the 
expression 

= dA& + ^ C [f ){i) AW A AW. (4) 

The integrability conditions 

dfW + ic ( ( ;» c) #AfW=0, (5) 

as well as the YM equations 

d*F(«) + Ic ( ^ (7) A^A*fW = 0, (6) 

are all satisfied. Here d is exterior derivative, A stands 
for wedge product and * represents Hodge duality. All 
these are in the usual exterior differential forms notation. 

Variation of the action with respect to the space-time 
metric g a b yields the EYM equations 



G a b — T ao . 
The stress-energy tensor is written as 

rp rpG i rpN 

1 ab — 1 ab + 1 abi 

where the gauge stress-energy tensor T^ b is 

(JV-l)(iV-2)/2 r 

ir a r b\ 2 

and the null dust energy-momentum tensor is 
T ab = *P(v,r)n a n b , 



rpG 

1 ab — 



E 



9ab 



(7) 
(8) 

,(9) 
(10) 



with ip{v, r), the nonzero energy-density and n a is a null 
vector such that n a = 5®, n a n a = 0. 

Expressed in terms of Eddington coordinates, the met- 
ric of general sphericall y sy mmetric space-time in N- 
dimensional space-times jl3l. [lH 25] is given by 



ds 2 



where 



(dSljv-s 



-A(v,r) 2 f(v,r) dv 2 +2eA(v,r) dv dr+r 2 {dn N - 2 ) 2 , 

(11) 



dej + sin 2 (6>i)d6^ + sin 2 (6»i) sm 2 (6 2 )d6 2 + 



' N-2 




\o 3 ) I del^ 



and {x a } = {v, r, 6\,... 6^-2}- For null dust, T vr 
must be nonzero and T£ = T£ for null energy condition. 
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This would imply A(v,r) — g(v) which could be set to 
1 without any loss of generality. It is useful to intro- 
duce a local mass function m(v, r) defined by f(v,r) = 
1 - 2m(v,r)/(N - 3)r( 7V " 3 ) (Hf. For m(v,r) = m(v) 
and A = 1, the metric reduces to the iV-dimensional 
Vaidya metric [15| . Therefore the entire family of so- 
lutions we are searching for is determined by a single 
function m(v,r) or f(v,r). 

It may be noted that in view of Eq. (3), the gauge 
field has only the angular components, Fg. s . with i =/= j, 

nonzero and they go as r~ 2 which in turn makes Ti? go 
as r~ 4 irrespective of N > 5. The null dust part will 
be given by T£ = ip(r,v). Note that the second term in 
Eq. (9) is evaluated to read 



Then the nonzero components would read as: T£ — 
i>{v,r), T% = T; = -(N - 3)(N - 2)Q 2 (v)/2r 4 and 

r£ = r£= ... = t*;- 2 = -(n-3)(n-6)q 2 (v)/2t\ 

It may be recalled that energy-momentum tensor (EMT) 
of a Type II fluid has a double null eigenvector, whereas 
an EMT of a Type I fluid has only one time-like eigen- 
vector m. 



(JV- 



l)(JV-2)/2 

E 



X(7 



(iV-3) (N-2)Q 2 (v) 



(12) 



For the energy-momentum tensor ((8|) and with the 
metric (TTTj) . the Einstein equations (O reduce to: 



il> = - 



(N - 2) df 



2r dv ' 
+ {N-3)Q 2 {v) =0, 



2 d 2 f , ,„ « 



»• ^ + (^-3)(2r^-(JV-4)(l-/) 



(N -3)(n-6)Q 2 (u) = 0. 



(13a) 
(13b) 

(13c) 



The last two equations are not independent and it suffices 
to integrate Eq. (|13b|) to give 



1 - 



1 - 



A/(-») 

r (N-3) 



nv- 



3)Q 2 (f) 
5 — 



(JV-5)r 



M{v) 2Q\v) ln(r) 



JV > 5, 



iV = 5. 



(14) 



where M(v) is an arbitrary function of v. Since YM T^ b 
go as r~ A (the same as for Maxwell field in N = 4), in- 
terestingly for all N > 5. That is why its contribution in 
/ will be the same for all N > 5 as in 4-dimensional 
Reissner-Nordstrom (RN) static or Bonnor- Vaidya ra- 
diating black- hole [26j . The nonradiating limit of this 
would be HD-Yaskin black-hole and not HD analogue of 
Reissner-Nordstrom. In other way round, this is making 
the Yaskin YM-black-hole radiate. 

There is however an important difference in the sign 
before Q 2 from the Maxwell case, which could be under- 
stood as follows [27|: gravitational potential $ at any r 
will go as 

$ = -(M-E(r))/r N - 3 , (15) 

where E(r) is the YM energy lying between r and infinity 
It is easy to compute E(r) 

(16) 



/"CO 

E{r) = / {Q 2 /r A )r n - 2 dr 

J r 



' (N — 5)r N ~ 5 



which would be negative for N > 5 and positive for N — 
4. This is what is responsible for the opposite sign for 
Q 2 as that for the Maxwell case in N = 4 (2||. So we 
have for N > 5 



M + (Q 2 r yv - 5 ) M Q 2 



„N-3 



(17) 



Unlike the Maxwell field, this is how YM field energy 
works in unison with mass and makes attractive contri- 
bution. For the Wu-Yang ansatz, the two fields are thus 
gravitationally distinguished. 

From Eq. (|13a|) . we obtain the energy density of the 
null dust with gauge charge as 



(AT-2) dM(v) , (N-2)(N-3) ^,, nA dQ(v) N>5 



ip(v,r) = 



2r 



W=2j—d7, 1 (AT_ 5 ) r 3 W(V)- 



dv 



2r 



(18) 



N = 5. 



and YM energy density and transverse stress are given 
by 



C(v,r) = (N-3)(N-2) 
P{v,r) = -(N-3)(N-6) 



Q 2 (v) 
2r 4 ' 

Q 2 (v) 
2r 4 ' 



(19) 
(20) 
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for TV > 5. Note that P(v,r) = for TV = 6 where 
it changes sign from positive to negative. The solu- 
tion obtained above represents a general class of non- 
static, TV-dimensional spherically symmetric solution of 
EYM theory describing a radiating black-hole. It con- 
tains TV-dimensional version of Vaidya [HI, [l9j] solution. 



In the static limit, (M = 0), it reduces to the black- hole 
solutions independently discovered in HD-EYM theory 




Also with (|14p. the Kretschmann scalar (JC = 
RabcdR abcd ) for the metric (HI]) reduces to 



K 



'^'yV w + 2A S:^? Q>)M{v) + 2{N -%%.^ +W) Q\v) TV>5, 



78 



M \v) + 24{12h p- 7) Q 2 (v)M(v) + 4(781n'(r)-841n(r)+31) Q4(7;)- 



TV 



(21) 



So the Kretschmann scalar diverges along r = 0. As ex- 
pected from the structure of stress tensor, the divergence 
in the last term is independent of TV which indicates that 
YM gauge field effect is dimension neutral. 



III. ENERGY CONDITIONS AND HORIZONS 



In the rest frame associated with the observer, 
energy-density of the matter will be given by, 



the 



i' = r v , c = -tj = -t; = 



(TV -3)(TV -2)Q 2 (v) 
2r 4 



(22) 

a) The weak and strong energy conditions (WEC/SEC): 
From Eqs. CED and (|2Ti]>. it is clear that > 0, C > 
always but P > only for TV < 6. That is the energy 
conditions will not in general be satisfied for TV > 6. For 
a Type II fluid, both WEC and SEC are identical, b) 
The dominant energy condition: It requires 



(23) 



which as before would be satisfied only for TV < 6. 

The line element of the radiating black-hole in 
Einstein- Yang-Mills theory has the form <jTTJ) with / 
given by Eq. (|14p and the energy-momentum tensor 
([8]). The luminosity due to loss of mass is given by 
Lm = —dM/dv, and due to charge by Lq = —dQ/dv, 
where Lm,Lq < 1. Both are measured in the region 
where d/dv is time-like [2!j[3(|. In order to further dis- 
cuss the physical nature of our solutions, we introduce 
the kinematical parameters. Following York [2{|, a null- 
vector decomposition of the metric is made of the form 



9ab 



-nJb - Lnb + Jab, 



(24) 



where. 



K,la=~f(v,rK + 6l (25a) 



W, (25b) 



lal a 



J[sin 2 (0 

n a n a = l a n a = -1, 
^7a6 = 0; 7a6 n h = 0, (25c) 



The optical behavior of null geodesies congruences is 
governed by the Raychaudhuri equation [29l . l30j 



dB 

dv 



KQ - R ab l 



ajb 



-e 2 



0~abO~ 



ab 



■ UJ ab UJ 



ab 



(26) 



with expansion 0, twist uj, shear cr, and surface gravity 
K. Here R a b is the TV-dimensional Ricci tensor, j c c is 
the trace of the projection tensor for null geodesies. The 
expansion of the null rays parameterized by v is given by 



6 = V a l a - K, 



(27) 



where the V is the covariant derivative and the surface 
gravity is 



K = -n a l b \7 b l a . 



(28) 



The apparent horizon (AH) is the outermost marginally 
trapped surface for the outgoing photons. The AH can be 
either null or spacelike, that is, it can "move" causally or 
acausally. As demonstrated by York [29| , horizons can be 
obtained by noting that (i) apparent horizons are defined 
as surface such that & ~ and (ii) event horizons are 
surfaces such that dQ/dv ~ 0. It follows that apparent 
horizons are the zeros of 



r (N-3) _ M „ (N_V Q 2 r (N-5) =0 7V>5; 



(29) 



M - 2Q 2 ln(r) = 



TV: 
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For Q 2 — 0, we have Schwarzschild horizon vah = 
(M) «-3 . In general Eq. admits a general solution 



r AH = t + -zr 



J^i? = exp 



Q 2 LambertWf- cxp( -ff 



/Q ) 



-M 



N = 6, 
N = 7, 

N = 5. 



(30) 



Here 



lent to the requirement that 



A = ( 4M ± 4a/M 2 - 4Q« 




FIG. 1: Plot of tah versus M and Q 



Thus for N = 5, the two roots of the Eq. (|29|) coincide, 
and there is only one horizon and for N > 5 there are 
two horizons, namely inner and outer horizons. Clearly, 
in the limit Q = 0, we have Schwarzschild horizon tab — 

On the other hand, the future event horizon is a null 
surface which is the locus of outgoing future-directed 
null geodesic rays, that never manage to reach arbitrarily 
large distances from the black-hole, and is determined via 
the Raychadhuri equation. It can be seen to be equiva- 



dv 2 



0. 



EH 



(31) 



The expression for the event horizon is the same as that 
for the apparent horizon with M and Q being respectively 
replaced by M* and Q* [II 13], where M* and Q* are 
effective mass and charge defined as follows: 

M*(v) = M(v)-^, (32) 



'(v) = Q(v)-^-. (33) 



IV. GRAVITATIONAL COLLAPSE 

In this section, we employ the above metric for investi- 
gation of YM gauge charge effect on formation of black- 
hole or naked singularity in collapse of null dust in the 
context of CCC. The physical situation is that of a ra- 
dial influx of gauge charged null dust in an initially empty 
region of the HD-Minkowski space-time. The first shell 
arrives at r — at time v = and the final at v = T . A 
central singularity of growing mass is developed at r = 0. 
For v < we have M(v) = Q(v) = 0, i.e., an empty 
HD-Minkowski metric, and for v > T, M(v) and Q 2 (v) 
are positive definite constants. The metric for v = to 
v = T is the HD-EYM Vaidya (discussed above), and 
for v > T we have the HD-EYM Schwarzschild [22|, 13] 
solution. 

In order to proceed further we would require the spe- 
cific forms of the functions M(v) and Q 2 (v), which we 
choose as follows Ujl [331 : 



and 



M(v) = 



Q 2 (v) 



0, v < 0, 

Xv(X > 0) < v < T, 
M (> 0) v > T. 



0, 



v < 0, 



fJ 2 v 2 (fi 2 > 0) < v < T, 
Q 2 (> 0) v > T. 



(34) 



(35) 
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Then the space-time is self-similar [3l|, admitting a ho- 
mothctic Killing vector 



d 



dr dv ' 
which is given by the Lie derivative 

L^ffab = £a;fc + 6>;a = ^9ab, 



(36) 



(37) 



where L denotes the Lie derivative. Let K a = dx a /dk be 
the tangent vector to the null geodesies, where k is an 
affinc parameter. Then 



g ab K a K b = 0. 
It follows that along null geodesies, we have 

CK a = rK r + vK v = C. 
Following [H}, we introduce 



p 

r 



and, from the null condition, we obtain 
M(v) (N-3)Q 2 (v) 



K" 



1 



r (JV-3) (N-5) r 2 



P 
2r' 



(38) 



(39) 



(40) 



(41) 



To study the singularity we employ the method devel- 
oped by Dwivedi and Joshi |32j. Consider the equation 
for future-directed outgoing radial null geodesies 



dr 
dv 



f 
2 



M(v) (N-3)Q 2 (v) 
r (N-5) r 2 



(42) 



The region / < is the trapped region and surface / = 
represent trapping horizon [32j. This is an ordinary 
differential equation with a singular point v = 0, r = 0. 
This singularity is (at least locally) naked if there are 
geodesies starting at it with a definite tangent. If no such 
geodesic exists, then singularity is not naked and strong 
CCC holds. To investigate the behavior near singular 
point, define 



y 



(43) 



Eq. ([32]). upon using Eqs. ([33]) and (|33"j). turns out 

to be 



dr 
dv 



-i \ JV-3 3) 2 2 

l - Xy -jN^f\ 



(44) 



If singularity is naked, there exists some value of A and 
fi, such that at least one positive finite value yo exists 
which solves the algebraic equation 



yo = lim y = lim — . 

1 — >0 v— >0 r— >0 v—*0 T 



(45) 



TABLE I: Variation of /j, c with D. For fj, > fj, c 
of collapse is a black-hole for all A (A > 10~ 12 ) 



the end state 



D = N 


Critical Value /j c 


Two equal Roots ~ yo 


6 


0.111111 


3.0 


7 


0.136081 


3.0 


8 


0.14907 


3.0 


9 


0.157133 


3.0 


10 


0.162648 


3.0 



Using (|44"|) and l'Hopital's rule we get 



y = lim y = lim — = lim — 

r— >0 u->0 r— >0 «->0 r r->0 V— >0 dr 



(46) 



1 ~ A vo yo 



which can be written in explicit form as, 



( iV - 3 ) ..2 7J 3 , X AN-2) 

(N-5 f Vo + Xy ° 



- yo + 2 = o. 



(47) 



This algebraic equation is the key equation which gov- 
erns the behavior of the tangent vector near the singular 
point. If the singularity is naked, Eq. (|47|) must have 
one or more positive roots yo, i-e., at least one outgoing 
geodesic that will terminate in the past at the singular- 
ity. While the absence of positive roots indicates that 
the collapse will always lead to a black- hole. Any solu- 
tion yo > of the Eq. (|4"T|) would correspond to the naked 
singularity of the space-time, i.e., to a future-directed null 
geodesic emanating from the singularity (v = 0,r = 0). 
The smallest such yo corresponds to the earliest ray ema- 
nating from the singularity and is called Cauchy horizon 
of the space-time. If yo is the smallest positive root of 
(|47l) . then there are no naked singularity in the region 
y < yo- It is easy to check that two roots of Eq. (|4T|) 
are always positive if A < X YM and the corresponding 
values of equal roots yo (for A = Ac) are shown in Ta- 
ble II. Thus, the occurrence of positive real roots implies 
that the strong CCC is violated, though not necessar- 
ily the weak CCC. The global nakedness of singularity 
can then be seen by making a junction onto HD-EYM 
Schwarschild space-time. 



In the absence of positive real roots, the central singu- 
larity is not naked (censored) because in that case there 
are no outgoing future-directed null geodesies from the 
singularity. It can be seen that collapse always leads to 
black-hole if [i > /i c (Table I) for all A < lfT 12 . In the 
limit fi — > (i.e. when gauge charge is switched off) 
our results reduce to those previously obtained by us for 
HD-Vaidya collapse [16], and in that case, Eq. (|47|) ad- 
mits positive roots for A < A c . Hence singularities are 
naked for A G (0, A c ], and censored (black holes) oth- 
erwise. Thus A = Ac is the critical value at which the 
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TABLE II: Variation of critical parameter A c and yo with D in HD-Vaidya and HD-EYM Collapse 



D = N 


HD-Vaidya Collapse 


HD-EM Collapse (/j, = 0.1) 




\v_ 1 
c N-2 


f N-3 ) 
y 2N-A j 




Double Roots(y = T$=r) 




Double Roots(yo) 


6 
7 
8 
9 
10 


27/2048=0.0132 
16/3125=0.00512 
3125/1492992= 0.00209 
729/823543= 0.00088 
823543/2147483648=0.00038 


8/3=2.6667 
5/2= 2.5 
12/5= 2.4 
7/3 =2.34 

16/7 = 2.28 


0.00238 
0.002135 
0.0009634 
0.00043546 
0.00018845 


2.91 
2.69 
2.56 
2.48 
2.42 



transition occurs, the end state of collapse switches from 
naked singularities to black holes. To conserve space, 
We have summarized the results in the Table I and II. 



All results concerning nakedness of singularity character- 
izing HD-Vaidya can be obtained from p — > (see [l6[ 
for details). 



J 



A. Strength of Naked Singularities: 

From the physical point of view, one of the most impor- 
tant features of curvature singularity is its gravitational 
strength. A singularity is termed gravitationally strong 
or simply strong, if it destroys by crushing or stretc hing 
any object which falls into it. A sufficient condition [331 ] 
for a strong singularity as defined by Tipler |34{ is that for 
at least one non spacelike geodesic with affine parameter 
k, in limiting approach to singularity, we must have 



lim k 2 tb = lim k 2 R ab K a K b > 

fc->0 fc->0 



(48) 



where R a b is the Ricci tensor. Eq. (j48|) . with the help of 
Eqs. <J34J) , g0]) and (gT]), can be expressed as 



lim k ib 

fe^O 



[N - 2)(JV-3) 



2(iV~ 5) 
x lim ((AT - 5)\y( N - 4) + Ap 2 y} 



kP 



(49) 



Our purpose here is to investigate the above condition 
along future-directed null geodesies coming out from the 
singularity. Eq. ([39j). because of Eqs. (J34]), <j35j), (gDJ) 
and (gj), yields 



P 



2C 



2- y + A^-2) + I^g M V 



(50) 



and geodesies are completely determined. Further, we 
note that 



1* = 

dk r r 



(51) 



which, on inserting the expressions for K r and K v , be- 
come 



dX 
~dk 



7JT= 2-y + Ay 



(JV 



_ 2)+ (iV-2) 2 *\P_ = C 



(iV-5) 



a* y 



2r 2 



(52) 



Using the fact that as singularity is approached, k — > 0, 
r — >• and X — > a + (a root of ([17])) and using l'Hopital's 
rule, we observe 



lim 

fc^O 



kP 



yo 

2 



(53) 



when limi 



>o f = Pq 7^ oo and hence Eq. (|49|) gives 
(A r -2)(A^-3) 



lim k ib = 



((^ - 5)Ay<; 



(AT-4) 



2(JV-5) 
4^ 2 y ) f > 0. 



(54) 



Thus along radial null geodesies strong curvature condi- 
tion is satisfied. Therefore, one may say that generically, 
the naked singularity is gravitationally strong [34J • Hav- 
ing seen that the naked singularity in our model is a 
strong curvature singularity, we also examine its scalar 
polynomial character. The singularity arising in the HD- 
Vaidya model was shown to be a strong curvature and 
also a scalar polynomial [l5|. The presence YM gauge 
charge does not affect this feature which is evident from 
the divergence of Kretschmann scalar. 



V. CONCLUDING REMARKS 

In conclusion, we have constructed nonstatic radiat- 
ing black-hole solutions of the coupled EYM equations 
for a null dust with gauge charge in HD, namely HD- 
EYM Vaidya. The HD-EYM Vaidya solutions are ob- 
tained by employing HD curved-space generalization of 
Wu-Yang ansatz [21] . Thus we have an explicit non- 
static radiating black- hole solutions of Einstein equations 
for non-Abelian gauge theory. This yields in 4D, the 
same results as one would expect for charge null dust 
in the Abelian theory (2(|, i.e., in 4D the geometry is 
precisely of the Bonnor- Vaidya form and the charge that 
determines the geometry is YM gauge charge. However 
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in HD-EYM radiating black-hole solutions deviate from 
Bonnor-Vaidya solutions because here the term Q 2 /r 2 in 
the solution (14) is dimension independent while it would 
go as Q 2 /r N ~ 2 for the latter. This is also reflected in the 
divergence of the Kretschmann scalar in Eq. (21). Note 
that the last term diverges as r~ 8 while for HD solu- 
tion with charge it would diverge as r^^ N ^ 2 \ This is 
the consequence of the Wu-Yang ansatz, which is also re- 
sponsible for charge working in unison (attractive) with 
mass in its gravitational contribution. 

The family of solutions discussed here belongs to Type 
II fluid. However, if M = Q = constant and the matter 
field degenerates to type I fluid, we can generate static 
black-hole solutions obtained in [22| by proper choice of 
these constants. In the static limit, this metric can be 
obtained from the metric in the usual spherically sym- 
metric form 

dr 2 

ds 2 = -f( r )dt 2 + — + r 2 (dn N ^ 2 ) 2 (55) 
f{r) 

by the coordinate transformation 

dv = A^y 1 { dt + e J^yj ■ ( 56 ) 

In case of spherical symmetry, even when f(r) is replaced 
by f{t,r), one can cast the metric in the form (TTTj) [35j |. 

We have also used this solution to study the end state 
of collapsing star and showed that there exists a regu- 
lar initial data which leads to naked singularity. The 



relevant question is what effect does the presence of 
the gauge charge have on formation or otherwise of a 
naked singularity. Our results imply that the presence 
of gauge charge leads to shrinking of the initial data 
space for naked singularity of the HD-Vaidya collapse. 
That is, it tends to favor black-hole. The gauge charge 
would contribute positively to gravity of the collapsing 
null dust. This should cover part of the parameter win- 
dow in the initial data set for naked singularity. This is 
what has been demonstrated. That is, the parameter set 
which gave rise to naked singularity in HD-Vaidya col- 
lapse may now lead to black-hole in the presence of the 
gauge charge. There exists a threshold value for fi, as 
shown in Table. II, the parameter window gets fully cov- 
ered ensuring formation of black-hole for all values of A. 
That is when fj, > fic the CCC is always respected. The 
important point is that collapse of null dust with gauge 
charge would favor black-hole in comparison to naked 
singularity. 

As final remarks it would be interesting to see how 
the results get modified in EYM theory with the Gauss- 
Bonnet combination of quadratic curvature terms and, 
in general, for the Lovelock polynomial |36| . 
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